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ABSTRACT 

The astrophysical site or sites responsible for the r-process of nucleosynthesis still remains an enigma. 
Since the rare earth region is formed in the latter stages of the r-process it provides a unique probe 
of the astrophysical conditions during which the r-process takes place. We use features of a successful 
rare earth region in the context of a high entropy r-process [S > lOOfc^) and discuss the types of as- 
trophysical conditions that produce abundance patterns that best match meteoritic and observational 
data. Despite uncertainties in nuclear physics input, this method effectively constrains astrophysical 
conditions. 

Subject headings: Nuclear reactions, nucleosynthesis, abundances 



1. INTRODUCTION 

Among many open problems in modern nuclear as- 
trophysics, a complete understanding of the rapid neu- 
tron capture or "r-process" remains one of the great- 
est challenges. The r-process is believed to be respon- 
sible for approximately half the heavy elements above 
the iron group. There are many promising candidate 
sites for this synthesis process, however the exact astro- 
physical site has not yet been identifie d; for reviews see 
( [Arnould et al. ||2007l, i Qian et ar]|2007[ ). Some examples 
ot candidate s ites includ' eTiigh-entropy supernovae ejecta 
(Meyer et ari992 Qian & Woosley|(l996| |Sumiyoslii ef 



aE 2 (JU(J[ I'hompson et al. 2001), compact object merg- 
ers (Latu mer et al . 1977; Meyer 1989; Freiburgh aus et 
al. 1999; Cori ely et al., 2005, Surman et al. 2008i1Mer 



zgcr et al. 2010 ), neutrino induced nucleosynthe sis in He 
sLells ( ,Woosley et al .ll990|IBaneriee et al.'2011'), or pos- 
sibly 7-ray burst scen arios (Surman & McLaughlin 2004j 



McLaughlin fc Surman||2005| jSurmaii et al.||2006| ). 



I'he study of the r-process is further complicated by 
the short lived nuclides far from stability. Recent experi 
ments have led to measurem ents on some nuclides partic 



ipating m t l 
Tones et al. 
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r-process e.g. 



Hosmer et al. (|2005| 



Baruah et al 
2010); Nishimura^t 



Mates 



Rahaman et al. 

(2009); Hosmer et al 

aTj ( |2011p ; |Quinn et al.| ( |201ip . However current exper- 
imentai dat a on neutron-rich isotop es is sparse; see e.g. 
figure 12 in Mumpower et al. (2011 ). Therefore, theoret- 
ical extrapola tions must be used as no experimental data 
are available ( Grawe et al.|2007 ). Much theoretical work 
is on going to understand the properties of nuclei far 
from stability and t h eir impact on the r-pr o cess, see e.g. 
Moller et al. (1995); Pearson et al. (1996); Rauscher et 
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(50X21. 
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(,2009,) 

( |2011| ); [Arcones &: Bertsch 



Pearson et a 



MumpOATOT 



Suzuki et al 



Suess & Urey ( 1956 ) demonstrated the existence of the 
= 180 and A 



195 peaks in an improved distribution 
of solar abundances. The prominent abundance features 
seen in this curve were quickly attributed cl osed neutron 
shells oc curring at magic neutron numbers (Burbidge et 

160) , lying away 
by a different 



al. 



19571). The rare earth peak {A 
Trom closed neutron shells, must occur 
mechanism. 

Three mechanisms have been proposed for the forma- 
tion of the rare earth peak. (1) It has been hypothesized 
that this peak may form by fission cycling if th e distribu- 
tion of fission fragments in (A, Z) is fa vorable ( Cameron 
" 19571 [Schramm et al.||197l| |Marti fc Zeh 1985; Goriely 



2011[ ). (2) [Surman et al.| ( I1997p discovered that the rare 
earth peak can form during treeze-out as material de- 
cays back to stability. Crucial to this argument was 
that this mechanism occur in a 'hot freeze-out', with 
temperature high enough to support significant photo- 
dissociation flow. This allowed material to 'funnel' into 
the peak in an attempt to sustain {n,j) <=^ (7,'t-) equi- 



librium. Otsuki et al. (2003) found similar abundance 
patterns provided the temperature was constant during 
freeze-out. (3) A third mechanism has also been found 
for lower tempe rature environments wh ere a 'cold freeze- 
out' may occur ( Mumpower et al.|2011 ). In this scenario 



the temperature dependent photo-dissociations flow is 
not present, however the rare earth peak may still form 
by a trapping mechanism involving slow neutron capture 
rates in the region. The formation process, in addition 
to being sensitive to freeze-out conditions was also shown 



al. 


2011 
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Arcones fc Martmez-Pinedo||201lf 



the most intriguing results to date regard- 
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ing r-process yields is that the stable elements observed 
in galactic halo stars (much older than the sun) with 
56 ^ 80, including the rare earth elements {Z = 57 
to Z = 71), are consistent wit h the scaled solar r-p rocess 
elemental abundances e.g. see'S neden et al.| (1996). This 
is fascinating because it suggests these elements were 
formed by an r-process mechanism which operates ove r 



a wide range of metallicity, e.g. see Cowan et al. (1997). 

An understanding of the production of the rare earth 
elements can be used to probe the freeze-out conditions 
of the r-process mechanism(s) which operates over a 
wide range of metallicity. This approach complements 
previous calculations which have focused on the condi- 
tions at the onset of the r -process i.e . the neutron-to- 



1997 



seed ratio (Hof fman et al. 1997 Meyer & Brown 
Freiburghaus et al. 199!^ 

In this paper, we introduce three new rare earth peak 
constraints on astrophysical conditions in the context of 
the site of the r-process. Using the formation of the 
rare earth peak as our key constraint, we isolate freeze- 
out conditions which are favorable for the production of 
these elements. We compare final abundances from sim- 
ulations to both the elemental abundances of halo stars 
and the isotopic abundances of r-process residuals. It is 
essential to note that (1) the r-process abundance pat- 
tern we observe may come about from a superposition 
of abundance patterns made from different conditions in 
the same astrophysical environment and (2) there exists 
significant variation between the final abundances using 
different nuclear models. We take this into considera- 
tion by averaging over and computing the variance of 
final abundances given different nuclear models and con- 
ditions. This allows us to compare final abundances, with 
error bars, to meteoritic and observational data. 

2. R-PROCESS MODEL 

Nuclides decay back to stability during freeze-out, the 
last stage of the r-process. To monitor this progress, it 
is useful to consider abundance weighted timescales for 
the primary reaction channels: neutron capture, photo- 
dissociation and /3-decay. 
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In equation [T] the neutron number density is denoted 
by Nn- The Maxwellian averaged neutron capture cross 
section for nuclide {Z,A) is {av)z.A and XjniZ,A) the 
photo-dissociation rate. The full /3-decay rate (including 
/3-delayed neutron emission channels) for nuclide [Z, A) 
is denoted by Xis{Z,A) and Y(Z,A) is the abundance of 
nuclei (Z.A). A reduced sum over the rare earth peak 
region (A = 150 to A= 180), denoted with a superscript 
"REP," may be taken when applicable. 

The neutron-to-seed ratio or K is another useful metric, 
we define as: 



R = 



(2) 



where F„ is the abundance of free neutrons. We con- 
sider 'low' neutron-to-seed ratio to be the time after neu- 
tron exhaustion, or i? = 1. 

The final abundances are highly dependent on the rate 
of decline in the temp erature and density during the last 
stage of the r-proc ess ( Arcones fc Martmez-Pinedo|20ir 



Mumpower et al.||2011" |. Among many possibilities, two 



thermodynamic evolutions are noteworthy and could oc 
cur in any astrophysical environment. A 'hot' freeze-out 
which proceeds under high temperatures (Tg > 0.6) when 
material decays back to stability and a 'cold' freeze-out 
which proceeds under low temperatur es (Tg 5-, 0.6) w hen 
material decays back to stability, e.g. Wanajo ( |2007 |. 

We contrast the differences in the evolution oi the rare 
earth region between the two thermodynamic trajecto- 
ries i n Figure [T] The traject ories used in this figure are 
from Mumpower et al. (2011 ) and the values in the time- 
line represent typical simulations which produce a satis- 
factory rare earth region compared to the solar isotopic 
abundances. We note that all of the simulations studied 
here shown qualitatively similar behavior. For further 



discussions of late time r-process dynamics see ( Surman 
let al."1997^'Surman fc EiIgeIl[200Tl |Arcones fc M 
IP- ■ " 



artmez- 



iPincdo 2011; Mumpowe r et aTp010[]20TT ) 

The rare earth region typically breaks from (rj,7) <=^ 
(7,n) equilibrium at Tg ^ 1, represented by point A in 
the figure. In the hot evolution, this is due to neutron 
exhaustion [R = 1) while in the cold evolution this is 
due to the rapidly dropping density. Once the neutron- 
to-seed-ratio is un ity, point B, the rare earth peak may 



begin to form (Surman et al. 
2011|). This process usually 



1997 Mumpower et al 
asts until p-dec&ys take 
over neutron captures in the region, point C. However, 
small changes in the pattern can continue after point 
C, until « a few T^-^^(not shown in the figure). 

The abundance averaged neutron separation energy is de- 
noted by (S'n). Following this quantity along the timeline 
we clearly see the rare earth peak forms as the nuclides 
decay back to stability (increasing values of (5'„)) and 
that cold freeze-out calculations typically extend farther 
from stability, closer to the neutron dripline than hot 
freeze-out calculations. 

In order to probe the conditions which are ideal for the 
r-process as matter decays back to stability, we study a 
range of possible outflow parameterizations. Our model 
consists of a monotonically decreasing temperature with 
density parameterized as: 



p{t) = piexp(-t/r) -I- p2 



(3) 



In equation [3j t denotes simulation time, pi -\- pi is the 
initial density at time t = 0, 3t = Tdyn with Tdyn the dy- 
namical timescale and A is a constant real number. For 
our calculations we take t ~ 27ms. The combination 
of r and entropy per baryon, 5, determines the distri- 
bution of seeds and the value of R at the onset of the 
neutron capture phase. We hold r constent and vary S 
in our calculations. The parameter n sets the thermo- 
dynamic behavior of the evolution at a low neutron-to- 



3 



seed ratio. A hot r-process evolution has typical values 
of n 2 which is character istic of wind models (Meyer 
2002| [Panov fc Janka ||2009[ ). Values of n > 5 are typical' 
of cola 7--process evol utions and correspond to a faster 
decline (Wanajo '2007"). 

Our r-process calculations are conducted with a nu- 
clear reaction network which includes approximately 
4000 nuclides relevant to the r-process. This code is 
fully dynamic in the sense that it does not assume any 
equilibrium conditions. The initial abundances for these 
nuclides are taken from an intermediate reaction network 



(Hix 



1999) with PARDISO solver (Schenk fc Gartner 



2004 1 . During the decay back to stability the primary re 
action channels are neutron capture, photo-dissociation, 
beta-decay, beta-delayed neutron emission, spontaneous 
fission and beta-delayed fission. 

Asymmetric fission is included during these freeze-out 
calculations ( jSeeger et al.,,1965i ) and operates for en- 
tropies above i>' = '200k b with proton to nucleon ratio 
or electron fraction. Ye = .30 and above 5* — 300^^ with 
Ye — -40. However, in these calculations, fission is not 
particularly important to the abundances at or above the 
atomic mass range of the rare earth peak until S 275fcs 
for Ye = .30 and S ~ 375fcs for Y^ = .40. 

To account for uncertainties in nuclear physics our 
nucleosynthesis calculations we use three different nu- 
clear data sets: F inite Range Droplet Mode l (FRDM) 



ih( 



Moller et al. 1995 



Rauscher et al . |2000 
1 Strutin sky Integra 



omas-i^er mi wit ^ ., 

(ETF SI-Q) ( [Pearson et al]|1996, ,Rauscher et al 



1998]), Extended 
and Quenching 



2000 



19981 and version 17 ot the Hartree F o ck Ijogoliubov 



moHe l (HFB-17) (Pearson et al.||2009 Goriely et al. 



2008 P We note that the /3-decay rates used for our 
calculations are only consistent with the FRDM nuclear 
data set. Since we arrive at our main conclusions after 
taking averages over and variances due to differences in 
nuclear models, we expect relatively small quantitative 
changes in our results from this approximation. How- 
ever, caution is warranted when attempting to judge the 
validity of individual nuclear models. The consistency of 
nuclear data for nuclear astrophysics investigations re- 
mains the subject of future research. 

3. COMPARING SIMULATIONS TO DATA 

In this section we compare final abundances to both 
halo star and solar data. We now comment on each of 
these data sets. 

Spectroscopic observations of galactic halo stars are 
useful in learning about the nature of nucleosynthesis 
processes that occur early in the history of the galaxy; see 
Sneden et al. (1996) and references within. The consis- 
tency ot the rare earth elemental abundance data among 
metal-poor halo stars is quite remarkable. This evidence 
suggests that these elements w ere s ynthesized in the same 
type of synthesis event ( jSneden et al.,,2008J and by a 
mechanism which operates over a wide ran ge of metallic- 
ity, for review see Cowan fc Sneden ( 2006 ) . Early galac- 
tic halo stars also provide an important diagnostic for 
r-process models as the s-process, occurring primarily in 
lower mass AGB stars, is not believed to have occurred 
yet. 

^ http://www.astro.ulb.ac.be/ 



We compare simulations to an averaged data set of 
five r-pro cess rich halo s tars: CS22892-05 2 (Sneden et 
[aDfl 996|, CS31082 -001 ( |Hill et al.|[2ro2t, HPT IMM 
(WestT xTet a'L||2000[), HD gmTO ( jiyans et al.||2006[ ), and 
BD+17"3248 ("cbwan et al.'2002).TTus averaged data set 
represents a typical r-proccss rich halo star. Before aver- 
aging the halo star data we arbitrarily normalized each 
halo star data to loge(Z = 63) = 0. We note that this 
choice of scaling is not critical to the anal ysis and has 
bee n used elsewhere in the literature, e.g. in |Roederer et| 



al.| (|2010p . We denote the averaged halo star elemental 
aoundances as H{Z). 

The solar isotopic abundance pattern corresponds to 
a superposition of many different events during the ori- 
gin of the solar sy stem (Kappel erl|1989[ Arlandini 1999 
|Cowan fc Sneden] [2006; Fodd ers et al |]2009] ). A t'ur- 



thcr complication arises when using solar data because 
in order to extract the solar r-process res iduals, s-process 



components mu s t be known accurately (Kappeler 1989 
Arlandini 1999 Roederer et al. 2010 1 . Despite these 



intricacies, solar r-process elemental yields for the rare 
earth region are strikingly similar to the elemental abun- 
dances of metal-poor, r-process rich halo stars. This 
means that solar data can also be used to constrain favor- 
able freeze-out conditions with the additional advantage 
of isotopic information. 
In the follow analy s is we use the most recent solar data 



from Fodders et al. (2009). This data set contains ac- 
curately measured abundances among the rare earth el- 
ements. To obtain r-process residuals, Nq^^j we take 
this solar data and subtract from it the stellar model 

The r- 



s-process abundances from Arlandini (19991. 
process residual s, Np,,^^ are in good agreement w ith pre- 
vious data sets ( |Kappeler]]T989l ]Arlandi"^|1999] ) . 



3.1. Definitions 



In this sub section (3.1 1 we introduce a series of defini- 
tions that we use in the remainder of section]3]to compare 
simulation to data. 

To compare simulations to halo star abundance data, 
H{Z), we first perform individual calculations using the 
three nuclear data sets (FRDM, ETFSI, HFB17) with 
conditions varying between < n < 10, lOOfc^ < S < 
400fcB, and with Yg = .30 or = .40. We then proceed 
to take an average of the final abundances from simula- 
tions over the different nuclear data holding the condi- 
tions constant (n,5',Ye). We then compare (1) element 
by element the rare earth peak {A = 159 to A — 167 
or Z — 64 through Z — 68) and (2) the ratio of total 
abundance in the rare earth peak to the platinum peak 
{A — 195 01 Z = 78) of the averaged simulation output 
to the corresponding values of our compiled halo star. 

To accomplish the first comparison we use G§^^ 
where, 



G 



REP _ 



H{Z)f 



E 



(4) 



In equation]! Y{Z) = Ea^(^'^) denotes the final 
abundance of me stable element with proton number Z 
and H{Z) is the average abundances of the five halo stars 
for given proton number. The summation is performed 
only over the rare earth peak elements as highlighted by 



4 



the superscript "REP" . A lower value of G^^^ implies a 
better match, with G§^^ = representing perfect agree- 
ment with H{Z). We scale the final abundance pattern, 

Y(Z), by a which is the solution to — — = 0. 

The A = 195 peak to rare earth peak ratio is an im- 
portant second constraint as both the lanthanides and 
heavier elements are believed to be produced in the same 
synthesis event and so their ratios must match observa- 
tional data. We compare the ratio of simulations to H{Z) 
using 772 where. 



aY{Z = 78) H{Z = 78) 



aY{Z) 



EzZllHiZ) 



(5) 



Only some of the halo stars have observations of plat- 
inum. So for the value of H{Z — 78) we average only 
those halo stars in which this element has been observed. 

To compare simulations to the solar r-process residu- 
als, NQ^ri-A), we proceed in a similar fashion using the 
previous two constraints. As an additional, third, con- 
straint we introduce and measure the amount of late time 
neutron capture occurring in the rare earth region as 
these nuclides decay back to stability. 

To gauge how close a simulation comes to reproducing 
the solar rare earth peak we use G^^^ , 



r<REP 



2^A=159 - 



(6) 



In equation IbI Y{A) = Y{Z, A) denotes the final 
abundance along an isobaric chain with atomic mass, A 
and NQ^r{A) is the solar r-process residual data also for 
an isobaric chain with atomic mass, A. Again, a lower 
value of G^^^ implies a better match, with G^^^ — 
representing perfect agreement with NQ_r{A). The final 
abundance pattern, Y(A), from simulation is scaled in- 



da 



^ 0. 



dividually by a which is the solution to 

Elemental abundance data from halo stars exhibit an 
A = 195 peak to rare earth peak ratio in agreement with 
the solar elemental abundance distribution. We thus use 
77^, the ratio of the two peak heights for solar data as 
well. 
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Our third constraint comes from the observation that 
too much neutron capture during freeze-out has negative 
conseque nces for the final r-process abundances ( Arcones 
[FMartmez-Pinedo|2011[|Mumpower et aLpoTTj ). I'his is 
particularly important tor the rare earth region, because 
the neutron capture rates in this region are on average 
faster tha n in other regions during the last stage of the 
r-process ([Arcones fc Martmez-Pinedo|2011 Mumpower 



et al.||201ip . Late-time neutron capture ettect shifts ma- 
terial troni below the peak {A — 150 to A = 158) to 
above the peak {A = 168 to A — 180) causing an over 
production of the heavier rare earth elements, which does 
not resemble either the solar r-process residuals, or halo 
star observations. 

We contrast the abundances of two simulations to ex- 
hibit the late-time neutron capture effect in Figure [2j In 



both panels of the figure we show a snapshot of abun- 
dances at different points during the calculation. The 
top panel shows abundances at i? = 1 (see point A, 
Figure [I]), the middle panel shows abundances during 
freeze-out (in between point B and C, Figure [l]) and 
the bottom panel shows the final isotopic abundances 
as compared to the solar abundances. The peak forms 
successfully in both cases, panel (a) with conditions 
(n — 6.0,5 — 350kB,Ye — .30) and panel (b) with con- 
ditions {n — 2.0,5 — 350fcB,Ye — .30). However, com- 
paring the middle and bottom snapshots in panel (a) we 
can clearly see the late time neutron capture effect dur- 
ing freeze-out. This effect is not present when comparing 
the same snapshots in panel (b). 

The late-time neutron capture effect is especially useful 
because, as we have just shown, it too is highly dependent 
on the conditions at late-times. We thus can use this ef- 
fect as a constraint on favorable r-process conditions. To 
estimate the amount of late time neutron capture which 
shifts material across the rare earth region we compute 
the difference between the value of Ga during freeze-out 
(between points B and C on the timeline) and the value 
of Ga using the final abundances. The summation now 
occurs over all rare earth elements {A = 150 to ^ = 180) 
signified by the superscript "REE" . 
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3.2. Determining Freeze-out Conditions 

In FigurelSlwe highlight astrophysical conditions where 
simulations oest match H{Z). The power law (n) from 
the density parameterization is plotted on the y-axis 
along with entropy per baryon (5) on the x-axis. Panel 
(a) simulations have Yg = .30 and panel (b) simulations 
have Ye = .40. Shown in each panel is a solid vertical line 
which denotes the point at which a sufficient neutron- 
to-seed ratio {R ^ 80 at Tg ~ 3) has been reached to 
prod uce a full r-process with elements out to A ^ 200 
( Hoff man et~aI1| 1997| [Meyer fc Brown ||1997| [Thompson 
|et al.| 2001 r We also denote the pomt at wnich hssion 
cycling influences the final abundances of the ^ > 150 
region by a dotted vertical line. 

In Figure[4[we highlight astrophysical conditions where 
simulations best match NQ_r{A). The axes and markings 
are the same as in Figure [3[ The left column shows re- 
gions of parameter space calculated with Y^ = .30 and 
the right column with Y^ = .40. 

We now discuss below each of the abundance features 
individually and their impact on constraining freeze-out 
conditions when comparing to halo star and solar data. 

3.2.1. Using Rare Earth Peak Formation 

Shown in both panels of Figure [3[ is a dark shaded re- 
gion (online red) where the scaled final abundances of 
rare earth peak elements from simulations best agree 
with the abundances of rare earth peak elements of 
H{Z). Since no single simulation reproduces exactly the 
H{Z) abundances, we must compute an appropriate cut- 
off for G^^^ . To this end, we measure the variance 
among the nuclear models by computing the standard 
deviation of the individual abundances at each value of 
Y(Z) given constant conditions {n,S,Ye). We addition- 
ally estimate the uncertainty in astrophysical environ- 
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riREP 



ment by computing a second average over the simula- 
tions. Propagation of error from averaging the nuclear 
models and averaging over astrophysical conditions pro- 
duces a value of G§ = 0.03. We thus take this value 
as a cut-off for shading the contour. 

Shown in both columns of panel (a) of Figure l4] is 
a dark shaded region (online red) where the scalea fi- 
nal abundances of rare earth peak elements from sim- 
ulations best agree with the abundances of rare earth 
peak elements of NQ^r{^)- We repeat the same proce- 
dure as mentioned above to compute the cut-off value of 
0.15 for shading this contour. 
Fission during freeze-out operates under higher en- 
tropy components {S > 200kB for = .30 and S > 
SOOfcs for Ye = .40) and influences abundances at or 
above the rare earth region at S ^ 275fcB for Yg — .30 
and S ~ 375fcs for Y^ — .40. We denote the point at 
which fission is influential in both Figures [3] and |4] by a 
dotted line in each case. Without fission, the dark shaded 
contour would end at or before the dotted line. As can be 
seen from the figures, the inclusion of fission extends the 
favorable calculated freeze-out conditions for rare earth 
peak formation. 

3.2.2. Using The A=195 Peak to Rare Earth Peak Ratio 

The light shaded region (online yellow) in both panels 
of Figure [3] represents the parameter space in which sim- 
ulations produce the same ratio of the A = 195 peak 
to rare earth peak as compared to the averaged halo 
star, H{Z). We take as an upper bound to this con- 
tour, rjz — 2, which represents an over production of 
the A = 195 peak relative to the rare earth peak by a 
factor of two, the value of 77^ « 2 comes from the ratio 

„"^£'f^I-^,'*^, . We note that the under production of the 

Z^Z=64 -Hl^) 

A = 195 peak is taken care of by the sufficient neutron- 
to-seed ratio constraint (solid vertical line). All together, 
this results in a constrai nt similar to the method used by 
(1997). Qualitatively similar behavior 
't and right columns of panel (b) of Fig- 



Meyer & Brown 
is shown in the le: 

ure [4] where calculated final abundances are compared 
to solar data. We take as an upper bound to this con- 
tour, r]A = 4: which again represents an over production 
of the A = 195 peak relative to the rare earth peak by 
approximately a factor of two. 

Interestingly, the left column of both figures, those 
with Ye — .30, show a second favorable region for this 
ratio at higher entropy. This is due to the oscillation of 
the A = 195 pe ak as a function of R due to fission cycling 



as discussed in Beun et al. ( 2006 
tween the two favorable regions t 



2008|). In the gap, be- 



195 peak is over 
produced by more than a factor of two compared to the 
rare earth peak. This effect also occurs in the Ye — -40 
case, however our simulations stop at 5* = 400fcB and so 
the effect is not shown. 

3.2.3. Using The Late-Time Neutron Capture Effect 

Returning to panel (c) of Figure |4j the shaded region 
(online green) in both columns represents the region of 
parameter space in which simulations do not exhibit a 
late-time neutron capture effect. We take as an upper 
bound to this contour, AG^^^ — 1, which represents 
significant fluctuations occurring across the rare earth 
region during freeze-out. 



Inspection of panel (c) reveals that the late-time neu- 
tron capture effect is sensitive to freeze-out conditions, 
this is particularly evident at high entropy in the left col- 
umn with Ye = .30. During a cold freeze-out (large values 
of n) the path moves far from stability where /3-delayed 
neutron probabilities are large. As the nuclides decay 
back to stability, /3-delayed neutron reactions provide an 
ample supply of neutrons for capture. Additionally, fis- 
sion actively introduces new seed nuclei into the region. 
Hence, the late-time neutron capture effect returns in 
this section of the flgure causing the shaded contour to 
disappear. At high entropy, in a hot freeze-out (lower val- 
ues of n), the late-time neutron capture effect does not 
occur (presence of the shaded region) due to the pres- 
ence of photo-dissociation and prolonged {n, 7) ^ (7, n) 
equilibrium. 

3.3. Further Analysis 

In order to emphasize the differences between the clas- 
sical neutron-to-seed ratio constraint (the region of pa- 
rameter space above the solid vertical line) and con- 
straints using the rare earth elements we highlight in 
each panel of Figures |3] and |4] a hatched region which 
represents our flnal constraint. 

In Figure [3] where simulation is compared with halo 
star data, the final constraint region is generally con- 
trolled by the A = 195 peak to rare-earth-peak ratio as 
the rare earth peak formation alone is a less stringent 
constraint. In Figure |4] where simulation is compared 
with solar data, the final constraint region is generally 
controlled by the same peak height ratio and addition- 
ally the late-time neutron capture effect. 

When comparing to both meteoritic and observational 
data given an initial Ye = .30, we find the entropy range 
favorable for rare earth peak formation is S" ^ 150A;b to 
S ^ 200fcB without fission. While the rare earth peak 
is formed by two different mechanisms in this entropy 
range, it does not favor one mechanism over the other. 
With the inclusion of fission a second favorable region 
appears around S ~ 325fcB to S ~ 375fcB in both fig- 
ures. In this region the rare earth peak formation and 
peak height constraints rule out the warmest and coldest 
scenarios when comparing to halo star data. When com- 
paring to halo star data at high entropy, hot conditions 
are disfavored due to the constraint on the rare earth 
peak elements and colder conditions are disfavored be- 
cause the A — 195 peak becomes over produced. When 
comparing to solar data, the late-time neutron capture 
effect rules out colder freeze-out scenarios (see discussion 
in previous section). 

We now explicitly investigate final abundances in the 
Ye — .30 case to the halo star and solar data. In order 
to observe the abundance pattern which is typical of the 
region suggested by the new constraint procedure we av- 
erage the final abundances in the constrained region and 
compute the standard deviation to measure the variance 
of these final abundances. 

The results of these calculations are shown in Figure 
[51 and |6] Panel (a) includes both averages over nuclear 
physics data sets and astrophysical conditions and the 
error bars represent variance from nuclear models and 
astrophysical conditions. The abundances in panel (a) 
correspond with simulations in the hatched region of the 
left column of Figure [3] and [4] respectively. Panels (b)-(d) 
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show the final abundances of individual nuclear models 
after averaging over their respective constraint regions 
and the error bars represent the variance from astrophys- 
ical conditions alone. 

As can be seen from Figure[5j these simulations do very 
well in reproducing the rare earth peak. When compar- 
ing to H{Z), the largest variation in the rare earth peak 
region comes from Terbium, Z = 65. Among all of the 
lanthanides, Praseodymium, Z — 59, shows the largest 
variance while Neodymium, Z — 60, is consistently un- 
der produced. Figure [6] also shows similar results for the 
rare earth peak elements. When comparing to the solar 
data, the largest variances in the final abundances can 
be found in the A ~ 180 to A 195 region. This ef- 
fect has been suggested to arise from a lack of long range 
corre lations in the nuclear models (Arcones & Bertsch 
20TT1) 



It is important to note in both figures that the variance 
seen in the final abundances of the rare earth elements is 
small. While subtleties persist in the individual nuclear 
models, global trends consistent with halo star and solar 
abundances distribution are present. As expected, we 
isolate conditions where the rare earth peak and the solar 
A = 195 to rare earth peak ratio mimic the solar data and 
further, prevent too much build up of the heavier rare 
earths. The success of this method is clearly dependent 
on how well we understand the input nuclear physics. 

4. SUMMARY AND CONCLUSIONS 

Using the delicate, out of (n,7) ^ {'^,n) equilibrium 
formation process of the rare earth peak we have stud- 
ied and isolated freeze-out conditions in high entropy 
{S > 100/cb) r-process calculations conducive for produc- 
tion of these elements. We compare the final abundances 
from simulations to both the elemental abundances of 
halo stars and the isotopic abundances of the solar r- 
process residuals. Both of these data sets suggest similar 
astrophysical conditions during freeze-out. For example, 
without fission, we find that favorable freeze-out condi- 
tions fall within a narrow entropy window {S ^ l^Oks to 
S ~ 200kB) given an initial electron fraction, = .30 
and timescale of t ~ 27ms. Given an initial electron 
fraction, — .40 and same initial timescale we find an 
entropy window of S ^ 300k b to S ^ 350fcs. 

The sensitivity of final abundances to uncertainties in 
nuclear physics is considered by including calculations 
with the FRDM, ETFSI and HFB17 nuclear models. 
We consider the variance among final abundances due 
to these three models and include it when determining 
the appropriate constraint for the formation of the rare 
earth peak (values of Gf^^ and G^^^). 

We use three abundance features to constrain astro- 
physical conditions: (1) We require the formation and 
correct shape of the rare earth peak {A — 157 to A = 167 
or Z — 64 to Z = 68). The importance of this process 
a nd its sensitivity to freez e-out conditions was discussed 



Mumpower et al. (2011). (2) We utilize the constraint 



on the ratio of the abundance in the rare earth peak 
to that of the A — 195 peak, as these elements are be- 
lieved to be synthesized in the same environment. (3) 
We utilize a new constraint on freeze-out behavior: the 
late-time neutron capture constraint. 

The late-time neutron capture constraint involves 
studying the shifting of material between isotopes dur- 



ing freeze-out in the rare earth region. We find the over 
production of the heaviest rare earth elements is highly 
dependent on the conditions during freeze-out, making 
it a valuable consideration for r-process models. This ef- 
fect can be measured by comparing abundances to solar 
r-process residuals as this data provides isotopic infor- 
mation. 

The impact of fission during freeze-out is to produce 
multiple islands of favorable r-process conditions, see e.g. 
the left column of Figures [3] or H] The late-time neutron 
capture constraint becomes tighter when fission occurs. 
This is due to two factors: (1) A cold r-process freeze- 
out extends far from stability where /3-delayed neutron 
probabilities are large, providing additional free neutrons 
to capture during freeze-out. (2) The inclusion of fission 
during freeze-out cycles additional material through the 
rare earth region which can capture the free neutrons. 
This suggests that if fission is important for the formation 
of the r-process elements, the astrophysical freeze-out 
conditions are more likely to be 'hot.' 

We have shown that the formation of the rare earth 
peak is a useful tool in determining favorable r-process 
freeze-out conditions. This diagnostic works well despite 
current uncertainties in nuclear physics. It complements 
the commonly used neutron-to-seed ratio and should be 
used in addition when analyzing astrophysical models. 
The size of the range of acceptable astrophysical condi- 
tions is dominated by uncertainty in the nuclear physics 
inputs. Reduction of this uncertainty through an im- 
proved understanding of the properties of nuclei far from 
stability and self-consistent nuclear data will increase the 
efficacy of this tool. 
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Figure 1. Contrasts the timeline of important events during tiie formation and evolution of the rare earth region between hot and cold 
evolutions. Point A represents the beginning of {rt,7) <^ (7i'^) freeze-out, point B neutron exhaustion (neutron-to-seed ratio of unity, 
R = 1), and point C when the timescale for /3-decays take over neutron capture. 




Figure 2. Highlights the late time neutron capture effect by showing the time evolution of rare earth region abundances versus the solar 
r-process curve (black) at three snapshots during the simulation. Panel (a) shows an r-process where late time neutron capture occurs 
using conditions {n = 6.0,5 = 350kg, Ye = .30). Panel (b) shows little to no late time neutron capture using conditions (n = 2.0,5 = 
350k B, Ye = .30). Calculations were performed with the FRDM nuclear data. Abundances in the top row are shown at i? = 1. Abundances 
in the middle row arc shown in between R = I and the point at which /3-decays take over. The bottom row shows final abundances. Refer 
back to Figure [l] 
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Figure 3. Shows how the rare earth region can be used to constrain astrophysical conditions when comparing simulations to typical 
halo star data. Calculations were performed with Ye = .30, panel (a) and Ye = .40, panel (b). The dark contour (online red) highlights 
conditions which produce a rare earth peak which best matches the average of five halo stars, H{Z). The light contour (online yellow) 
highlights conditions which produce an A = 195 peak to rare earth peak which best matches H{Z). The intersection of these constraints 
is showing by the hatching. A full r-process out to A ~ 200 occurs above entropies denoted by the solid vertical line. Above the dotted 
line fission influences the abundances of A > 150 nuclides. 
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Figure 4. Shows how the rare earth peak can be used to constrain astrophysical conditions when comparing simulations to the solar 
r-process data. The axes are the same as Fig. [3] The left column displays computations performed with Ye = .30 and right column with 
Ye = .40. Panel (a) highlights regions (online red) where the rare earth peak elements A = 159 to yl = 167 from simulations best agree 
with the rare earth peak elements of the solar r-process data. Panel (b) shows regions (online yellow) where simulations produce the same 
ratio of the A = 195 peak to rare earth peak as compared to the solar data. Panel (c) highlights regions (online green) where the late time 
neutron capture effect is minimal. The intersection (shaded region) of these three constraints are highlighted in each panel and repeated 
in Panel (d) with all three constraints. A full r-process out to A ^ 200 occurs above entropies denoted by the solid vertical line. Above 
the dotted line fission influences the abundances of yl > 150 nuclides. 
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Figure 5. Shows calculated final averaged elemental abundance patterns (solid red) along with uncertainties compared to the average 
of five halo stars, H{Z) (dots). All calculations were performed with Ye = .30. Panel (a) shows an average over 3 nuclear data sets 
and astrophysical conditions from the constraint region of Figure [s] Panels (b)-(d) show the average over astrophysical conditions in the 
constraint region for each individual nuclear model. 
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Figure 6. Shows calculated final averaged isotopic abundance patterns (solid red) along with uncertainties compared to the solar r-process 
residuals (black dots). All calculations were performed with Ye = .30. Panel (a) shows the average over 3 nuclear data sets and astrophysical 
conditions from the constraint region of Figure |4] Panels (b)-(d) show the average over astrophysical conditions in the constraint region 
for each individual nuclear model. 



